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In the framework of light-cone QCD sum rules, we study the valence quark distri- 
bution function q(xg) of a pion for moderate xg. The sum rule with the leading 
twist-2 wave function gives q(xg) = Vtt(^s)- Twist-4 wave functions give about 
30% for xg ~ 0.5. It is shown that QCD sum rule predictions, with the asymptotic 
pion wave function, are in good agreement with experimental data. We found that 
a two-hump profile for the twist-2 wave function leads to a valence quark disribu- 
tion function that contradicts experimental data. 



1 Introduction 

Two types of objects in hadron physics are used to represent quark and gluon 
distributions in a hadron. One is the quark and gluon distribution functions 
that appear in deep-inelastic lepton-hadron scattering, and the other is the 
light-cone wave functions of hadrons introduced in pQCD by Chernyak and 
ZhitnitskyB to describe hadron form factors at large Q 2 . In this paper, we 
investigate a relationship between these objects. To make a prediction for the 
valence-quark distribution in terms of light-cone wave functions, we combine 
the ideas of calculating the deep-inelastic_scattering amplitudes in the QCD 
sum rule approachH as suggested by IoffeS with the so-called light-cone QCD 
sum rulesu. 

In this work we propose a method to calculate the pion structure func- 
tion directly in terms of light-cone wave functions. The starting point of this 
method was formulated in Rcfa, where it was pointed that if xb = Q 2 /(2pq) 
is not close to the boundary values x = and x — 1, then the imaginary part 
of the deep-inelastic scattering amplitude is determined by small distances in 
the i-channel. A general proof of this statement follows from the fact that at 
large \p 2 \ and |g 2 |, the nearest singularity in the i-channel is at t = — l x _ B J . B V 1 
in the kinematics that is used in this paper. Thus, in the case of intermediate 
i^we can use the OPE to construct QCD sum rules. 
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2 Correlator 

To calculate the quark distribution function for a pion, we consider the follow- 
ing correlator: 

T ppX (p,q,k) =-ij d 4 xd 4 ze^ x+ ^ < 0\T{]l(x),j d p (z), 3 d x (0)\n~(k) >, (1) 

where A; is a pion momentum, and where 

il = uwnd, j* = dj p d, (2) 

k 2 = 0, q 2 = (p + q-k) 2 , t=(p-k) 2 = 0, S = (p + q) 2 , Q 2 = -q 2 . (3) 
^From relations (||) it is easy to determine that 

{2k,p + q) = s + Q 2 , (2pk)=p 2 . (4) 
We calculate the discontinuity in s at fixed p 2 and Q 2 of the correlator ([!]), 

ImT^x = i [T ppX (p 2 , q 2 , s + is) - T ppX (p 2 , q 2 , a ~ ie)] , (5) 

where p 2 and q 2 are space-like vectors, p 2 < 0, q 2 < 0, such that \p 2 \, \q 2 \ 3> 
Aqcd- In the scaling limit, we assume that \p 2 \ <C \q 2 \ and keep only the 
first nonvanishing terms in an expansion in powers of p 2 /q 2 . Perturbative 
logarithmic corrections are taken into account by the Alterelli-Parisi equation. 

We calculate ImT ppX in the physical region of the s-channel, and the pion 
contribution in this amplitude has the following form: 

ImT p . pX =pJ-^Imlj d 4 ze^ z < 7r(p)\T{j d p (z), j d (0)}\n(k) > j . (6) 

The general form for Im{J d i ze lqz < ir(p)\T{j d (z), j d (0)}\Tr(k) >} is 

Ai{s,Q 2 )(q 2 g p \ - q p 2) q x } ) 



+A 2 ( S ,Q 2 )(q 2 9pX -q^ki 1) ~ qp 2 k x 2) +q ( p 1) q x 2) ) 
+ B l{ s,Q 2 )( P - P -C q u) 



r / p \ q 2 



q 
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where q^ 1 ) = g and g( 2 ) = p + q — k are the momenta of the virtual photons; 

(g(l))2 = (g(2))2 = g 2 j t ig dear tha< . 

4x 2 B q d (x B ) = ^(B 1 ( S ,Q 2 ) + B 2 { S ,Q 2 ) + B 3 ( S ,Q 2 ) + B i ( S ,Q 2 )) Q 2 -> oo (8) 
where 

/m|y d 4 ze^ < 7r(p)|r{ Jp d (z),jt(0)}|7r(p) > 
x 2 B q d (x B ) ( pq \ ( pq \ 

Here q d (xB) is (i-quark distribution function of a pion. 

To find the quark distribution function in this paper, wwe consider the 
tensor structure p^p p p\ in correlator (Q). We define the imaginary part of the 

correlation function for these tensor structures as ^ B t(p 2 , x B ). We can write 
the following dispersion relation for this function: 

t(p 2 ,x B) =u(^ + j ds ^ ds y do) 

where the second term in the right-hand side is the higher-states contribution. 
To supress the exited states contribution, as usually done in QCD sum rules, 
we will consider instead of t(p 2 ,x) its Borel transform in p 2 : 

t(M 2 , XB ) = { -^f- (^yv,^) 

= ~U {x 2 B q d {x B ) + J p( S , XB )e-^ M2 d S ^j . (11) 



3 QCD sum rule 

To find q d (x B ), we apply the OPE near the light-cone x 2 = 0. In contrast to 
conventional QCD sum rules based on the OPE of a T-product of currents at 
small distances, we apply an expansion near the light-cone expressed in terms 
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of nonlocal operators, i.e. matrix elements that define hadron light-cone wave 
functions of increasing twist. The amplitude of two-quark nonlocal operator 
has the following form: 



< 0\u(x) W d(0)\7r(k) >= [ due^ kxu ( V7r {u) + x 2 9l {u) + 0(z 4 )) 

Jo 

+U (z M - ^k^j J due- lkxu (g 2 (u) + 0(x 2 )), (12) 

where ip n (u), and g±(u) and g 2 (u), are the leading twist-2 and the twist-4 pion 
light-cone wave functions, respectively. 

The leading twist-2 operator gives the following contribution to t(p 2 ,x) 

(13) 

Using eq.(|lO|) for t(p 2 ,x B ), the leading twist-2 wave function contribution is 

q d {x B ) = <Pk( x b)- (14) 

It is clear that this relation corresponds to a pure parton picture in which the 
pion consists of two valence quarks only: 

l ,i 

q d {x B )dx B = 1, / x B q d (x B )dxB = 1/2. (15) 
Jo 

These relations follow from the normalization of the twist-2 pion wave function 
( L 1 tp % (v,)du =1), and from its symmetry: ^(u) — ip T (l — u). 
In this paper we use the asymptotic pion wave function: 

<p v {u) = tp {asym: \u) = 6u(l -it). (16) 

Note that there are alternative models for ip w : the Chernyak-Zhitnitsky wave 
functions ip^ z and the function that was introduced in RefQ. Because of their 
two-humped shape, these alternative models for the pion wave function lead 
to poor agreement between the QCD sum rule prediction and experiment. 

To make a reasonable prediction, we have to estimate the twist-4 wave 
function contribution. In this paper, we take into account the contribution of 
the two-particle wave functions of twist-4. This contribution to t(p 2 , x B ) is 

U ft s 4/V ( gi(x B ) + G 2 (x B ) 1 d gi {x B ) \ 
V^ XB) = ^ { T B + 2 92{XB) - -iLxV ) ■ 
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where G 2 (u) = - Cdvg 2 (v). 

Using eqs . ( |To| , [l3| , [T6| ) , after making a Borel transformation, we have 

/ \ 4 (gi{x B ) + G 2 {x B ) , 1 , . dg 1 (x B )\ 

= q d (x B ) + C{x B )e- m ^' M2 . (18) 

Here, we take into account the contribution of the Ai-meson to evaluate the 
higher states contribution. 

To estimate twist-4_contributipn, we use the following set of twist-4 pion 
wave functions (see RefB and Refill): 

5 1 

gi(u) = S 2 - + -eS 2 (uu(2 + I3uu 
+2u 3 (10 - 15u + 6u 2 ) ln(u) + 2m 3 (10 - 1M + 6m 2 ) hx(u)) , 

92{u) = ^-8 2 uu(u - u), (19) 

where u = 1 — u. One of the parameters is defined by the matrix element 

< Q\g s uG af3 ~fi3u\ir(p) >= i5 2 fnP^- (20) 

The QCD sum rule estimate of Refi yelds S 2 = 0.2GeV 2 at fi ~ IGeV. At 
H ~ IGeV, e ~ 0.5 (see RefQ). 

To fix the function C{x B ), we take the QCD sum rule ( |l7| ) in the limit 
M 2 — > oo to be valid up to terms \/M 2 . This condition leads to the relation: 
C(x B ) = &gal. 

Finally we obtain the following QCD sum rule: 

q d (x B ) = <p«{x B ) - U{x B ) (jL + -L-e- m M/ M ^ . (21) 

For M 2 ~ IGeV 2 and 0.2 < x B < 0.8, the contribution of A i -meson, 
which imitates the higher state contribution, and the contribution of twist-4 
wave functions, is less than 30%. The results are depicted in Fig.l. 

Note that in our calculations we do not take into consideration the contri- 
bution of leading logarithms. So, in order to make a comparison with exper- 
iment, we have to use the quark distribution function at low Q 2 where these 
logarithms are small. Evolution of experimental data to low Q 2 has been done 
in Refu, and we use the results of this paper as experimental data. 

It is interesting to estimate the second momemnt of the quark distribution 
function. Assuming that the region near the end points x B = 0, 1 (where our 
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Figure 1: 



consideration is not valid) gives a small contribution to the second moment, 
the sum rule (|2l"| ) gives the following sum rule for the second moment: 



The analysis of this sum rule gives — 0.27 ± 0.05 at low normalization 
point, which is in a good agreement with experimental data: M% — 0.3 (see 



4 Conclusions 

The light-cone QCD sum rules considered in this paper lead to satisfactory 
agreement with experimental data. To improve the prediction of the QCD 
sum rule, it is nessesary to take into consideration twist-4 quark-gluon wave 
functions. One should note that the present experimental data for the quark 
distribution functions was obtained from an analysis of the Drell-Yan process, 
which is uncertain through the so-called K-factor. The choice ip n — cp( as J' m ) 
gives the best agreement between the QCD sum rule prediction and the ex- 
perimental data. Alternative models for the pion wave function having a two- 
humped profile lead to a poor agreement with the experimental data. 
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